Enhanced diffraction by a rectangular grating made of a negative 
phase-velocity (or negative index) material 
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Abstract. The diffraction of electromagnetic plane waves by a rectangular grating formed 
by discrete steps in the interface of a homogeneous, isotropic, linear, negative phase- 
velocity (negative index) material with free space is studied using the semi-analytic C 
method. When a nonspecular diffracted order is of the propagating type, coupling to that 
order is significantly larger for a negative index material than for conventional material. 
The computed coupling strengths reported here are in agreement with recent experiments, 
and illustrate the role of evanescent fields localized at the grating interface in producing 
this enhanced coupling. 
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The study of electromagnetic fields received a major boost in 2001 with the experimental 
confirmation of negative refraction by dielectric-magnetic materials with negative real 
permittivity and negative real permeability When considered as linear, isotropic and 
homogeneous, these materials possess a complex-valued refractive index whose real part 
is negative, and are therefore often called negative index (NI) materials. Alternatively, 
these materials can be referred to as negative phase-velocity (NPV) materials, as the phase 
velocity and the time-averaged Poynting vectors therein are antiparallel. Other names are 
also in circulation; but, regardless of the name used, the common observable phenomenon 
is negative refraction ^1 . 

All experimental realizations of NPV materials thus far are as periodically patterned 
composite materials. The unit cell comprises various arrangements of conducting filaments 
in order to realize both dielectric and magnetic response properties in the same frequency 
range. For such a material to be considered as effectively homogeneous, the unit cell must 
be electrically small — i.e., it must be considerably smaller than the free-space wavelength 
as well as the wavelength in the material In the wedge sample used by Shelby et al. 
PP, the unit cell size was about one-sixth of the free-space wavelength. The finite size of 
the unit cell meant that one of the two exposed surfaces of the sample was not planar, but 
rather formed a rectangular grating jSJ. Because of this grating, the specular (or zeroth- 
order) reflected/refracted plane wave could be accompanied by nonspecular (higher-order) 
diffracted plane waves [Sj. 

Recently, using Ansoft's HFSS package, an electromagnetics equation solver based on 
the finite-element method. Smith et al. [7] examined the plane-wave response of a wedge 
having the same properties as that used by Shelby et al. One of the two exposed 
surfaces was set up as a shallow rectangular grating. The simulation revealed not only 
a zeroth-order transmitted plane wave but also, unexpectedly, a strong first-order trans- 
mitted plane wave. The latter was much stronger for a NPV wedge than for a wedge with 
identical dimensions but made of a positive phase-velocity (PPV) material. Although 
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higher-order refracted plane waves were not experimentally observed by Shelby et al. 
the first-order was definitively present in later experiments jlj [7] . 

While full-wave simulations confirm the phenomenon of enhanced diffraction at the 
NPV grating, the numerical results are not of use in determining the physical origin of the 
effect. An alternative approach is furnished by mathematical treatments of diffraction by 
gratings based on analytic field expansions. These treatments have continued to develop 
over the last hundred years, and have now acquired a considerable degree of sophistication. 
We apply one of these methods here, the so-called C method, both to independently verify 
the simulation results obtained by the finite-element method and to provide further insight 
as to the underlying mechanism 7 . The C method, originally developed for dielectric 
gratings [SI, was modified to handle dielectric-magnetic gratings 0. 

In a rectangular coordinate system {x,y,z), we consider the periodically corrugated 
boundary y = g{x) = g{x + d) between vacuum and a homogeneous, isotropic, dielectric- 
magnetic material, with d being the corrugation period, as shown in Figure ^ The region 
y > g{x) is vacuous, whereas the medium occupying the region y < g{x) is characterized 
by complex-valued scalars £2 = e2R + ie2i and iJ,2 = fJ^2R + 11^-21, such that £2/ > and 
fJ'21 ^ 0. The refracting medium is assumed to be effectively homogeneous at the angular 
frequency of interest. A linearly polarized electromagnetic plane wave is incident on this 
boundary from the region y > g{x) at an angle ^Oi (I^ol < "^Z^), with respect to the y axis. 

Let the function f{x,y) represent the z-directed component of the total electric field 
for the s-polarization case, and the z-directed component of the total magnetic field for 
the p-polarization case ^0]. Outside the corrugations, f{x, y) is rigorously represented by 
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Rayleigh expansions as 



f{x, y) = exp \i {uqx - P^^^y) 
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Here, 
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are scalar coefficients to be determined; and 
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where c is the speed of hght in vacuum and u) is the angular frequency. Note that /3n^ is 



either purely real or purely imaginary; and the conditions Re 
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> and Im 



> 



Vn are appropriate for plane waves in the vacuous half-space y > max g{x). The refracted 



plane waves must attenuate as y — > —00, imposing the condition Im 
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> 0. Fulfillment 



of this condition automatically fixes the sign of Re 
and fi2R- 



, regardless of the signs of e2R 



After implementing the C method on a computer, the coefficients {/5n}„=t^ and 
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are determined. Diffraction efficiencies en = (Re (3n^ /Po^^) are defined 



for the reflected orders. If dissipation in the refracting medium is small enough to be 



ignored, diffraction efficiencies e'^ = (Re 



(2) 



/a (3q^^) [r^l^ are defined for the refracted 



orders, where either a = fi2 (s-polarization) ov a = €2 (p-polarization). 

The grating of Shelby et al. (IJ has a rectangular profile, shown in Figure^ with the long 
and the short sides in the ratio 3 : 1 and period d = 15.81 mm. To simulate this profile for 
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the C method, g{x) was replaced by the truncated Fourier sum "^^=0 sin{2Trnx / d+ ipn) , 
which was found adequate to represent the grating profile. 

In Figure[21 we show the diffraction efficiencies obtained as functions of a frequency 
of 11.75 GHz, the highest frequency used by Smith et al. 7'., which corresponds to 
ujd/c = 27r/1.58. The refracting medium is of either the the NPV (e2 = —5 + iO.Ol, fj,2 = 
— 1 + iO.Ol) or the PPV (e2 = 5 + i0.01, /U2 = 1 +«0.01) type. Calculations were made for 
both the s~ and the ^-polarization incidence conditions. As confirmed by perturbation 
analysis jl2j . in this case the transformation {e — e*, fi — > —fJ.*} does not greatly affect 
Cq (Fig. I^K); except at low j^ol- In contrast. Figs. [5)3 and|2t show that the nonspecular 
reflection efficiencies e^^, are significantly affected by the type of the refracting medium. 
In particular, the interplay between the polarization state and the angle of incidence 
leads to (i) very little difference for the PPV grating between the two polarizations; (ii) 
considerable difference for the NPV grating between the two polarizations, and (iii) a shift 
of the difference between PPV and NPV gratings from one polarization to the other as the 
sign of the angle of incidence changes. As the refracting medium is dissipative, diffraction 
efficiencies for the refracted orders cannot be defined. 

To study refraction efficiencies explicitly, we made the refracting material nondissipative 
(e2 = ±5, fi2 = il)) while the other parameters were kept the same. The higher-order 
reflection efficiency curves display the same qualitatitive behaviors as in Figure |2J and are 
therefore not shown in Figure 13 All refraction efficiencies are greatly affected by the type 
of the refracting medium, as becomes clear from Figs. OK^c. In particular, the coupling 
of power into the n = ±1 refracted orders (Figs. |3h s^nd c) for the NPV grating is much 
larger than for a PPV grating, a fact that again is in agreement with the results of Smith 
et al. (7|. 

A similar conclusion is drawn from Figure^ for which the efficiency curves were com- 
puted at a frequency of 9.0 GHz (the smallest frequency used by Smith et al. [7]), the 
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other parameters being the same as in Figure |31 Whereas the zeroth-order (Fig. the 
only propagating order reflected at this frequency) reflection efficiency is much less af- 
fected by the type of the refracting medium than in our previous examples, the refraction 
efficiencies depend strongly on the refracting medium being of the PPV or the NPV type. 

The increased coupling into higher-order diffracted plane waves represents an important 
distinction between the behaviors of PPV and NPV materials, and confirms that surface 
periodicity plays a much more significant role in the latter [7j. The physical origin of 
the enhanced higher-order reflection efficiencies for the NPV grating can be understood 
in an approximate manner by examining the reflection efficiencies for shallow gratings. 
In general, a grating couples the specular coefficients (po fmd tq) with all nonspecular 
coefficients (/)„ and t„, |n| > 0) ;6], thereby coupling Bq with all en, \n\ > 0. All coefhcients 
for a grating of arbitrary profile are obtained together by some numerical technique. 

In the limit of the grating profile becoming planar, the reflection coefficients (and, 
therefore, the efficiencies) coincide with those for a planar interface. The introduction of 
shallow corrugations only slightly perturbs the coefficients |13j . and we may then write 

(2) 

In the limit of small dissipation, f3n is real-valued for propagating orders in the refracting 

(2) 

medium. As this medium transforms from PPV to NPV, (3n changes sign along with a, 

(2) 

and the magnitude of pn is thus unchanged. But for evanescent (nonspecular) orders, 
is imaginary, and thus always positive in accordance with causality. Hence, for nonspecular 
orders, the reflection coefficients are inverted as the refracting medium is transformed from 
PPV to NPV. The evanescent orders then, which play a minor role for the PPV grating, 
mediate a much larger interaction between the specular and the nonspecular modes for 
the NPV grating. 

Equation (jlj) shows that there is the possibility of a pole for one or more reflection 
orders. A pole indicates the presence of a bound surface mode, identical to the surface 



plasmon that occurs at the interface with a metal at optical wavelengths. At least in 
a perturbative sense, the enhanced diffraction by a NPV grating has the same origin 
as other plasmon-mediated effects, including enhanced light transmission ^1] and the 
"perfect lens" effect [T^ . 

The detailed treatment presented here of diffraction from a NPV grating reveals en- 
hanced diffraction, in agreement with recent experiments. As NPV materials are consid- 
ered for various applications, the results found here indicate that great care is necessary 
in describing the surfaces of NPV materials. In particular, the recent artificially struc- 
tured metamaterials are based on periodic cells, whose inherent periodicity can lead to a 
quite significant surface nonhomogeneity, which may assume even greater importance for 
surfaces that are not nominally planar. 
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Figure 1: Schematic of the diffraction problem solved. The refracted plane waves are shown 
as if the refracting medium is of the NPV type. The specular reflected and refracted orders 
are denoted by n = 0, while nonspecular orders are denoted by n 7^ O.The inset shows the 
shape of the grating profile used by Smith et al. [7j as well as for Figures I^HS 
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Figure 2: Computed diffraction efficiencies as functions of tlie angle of incidence when 
u)d/c = 27r/1.58, for both p- and s-polarized plane waves. The refracting medium is either 
of the PPV (£2 = 5 + iO.Ol, /U2 = 1 + iO.Ol) or the NPV (e2 = -5 + iO.Ol, H2 = -1 + iO.Ol) 
type, (a) eg; (b) e^^^; (c) e^J.^. 
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Figure 3: Computed diffraction efficiencies as functions of the angle of incidence when 
u)d/c = 27r/1.58, for both p- and s-polarized plane waves. The refracting medium is either 
of the PPV (e2 = 5, /X2 = 1) or the NPV (e2 = -5, /X2 = -1) type, (a) e^; (b) eL^; (c) 
el,. 
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Figure 4: Computed diffraction efficiencies as functions of tlie angle of incidence 6q, when 
ujd/c = 27r/2.087, for both p- and s-polarized plane waves. The refracting medium is 
either of the PPV (e2 = 5, H2 = 1) or the NPV (e2 = -5, fi2 = -1) type, (a) eg; (b) ej; 
(c) eLi; (d) el,. 
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